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Abstract 

We analyze subtracted dispersion relations for meson correlators at finite baryon 
density and temperature. Such relations are needed for QCD sum rules. We point out 
the importance of scattering terms, as well as finite, well-defined subtraction con- 
stants. Both are necessary for consistency, in particular for the equality of the longi- 
tudinal and transverse correlators in the limit of the vanishing three-momentum of 
mesons relative to the medium. We present detailed calculations in various mesonic 
channels for the case of the Fermi gas of nucleons. 
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1 Introduction 



There has been a continued interest in the studies of in-medium properties of 
hadrons. As the QCD ground state is largely modified at finite temperature 
and/or baryon density, one expects that masses, widths, coupling constants, 
and other characteristics of hadrons change considerably [1,2]. In particular, 
the properties of light vector mesons are of the special interest, since their 
modification may influence dilepton spectra measured in relativistic heavy-ion 
collisions. In fact, the observed excess of dilepton pairs in such experiments 
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[3,4] is commonly explained by the assumption that the masses of the vector 
mesons decrease or their widths become larger, or both [5-8]. 

Different theoretical methods have been used to analyze medium effects on 
hadron properties. Let us mention effective hadronic models [9,10], low-density 
theorems [11,12], quark models [13], and QCD sum rules [14-25]. The results 
of this paper are relevant to this last method. In QCD sum rules one needs 
the in-medium dispersion relation for the meson correlator. Furtheremore, this 
dispersion relation is subtracted in order to improve the reliability of the sum 
rule method. We point out that the existent calculations [16-18,21-24] have 
not been very careful in this issue. They include the contributions form the 
so-called scattering term in the spectral density, but overlook the fact that in 
addition there is a well-defined, finite subtraction constant in the dispersion 
relation that has to be kept. This subtraction constant, which depends on the 
channel, is crucial for consistency of the approach. For instance, only when it 
is kept then the transverse and longitudinal correlators in vector and tensor 
channels are equal in the case of the vanishing three-momentum in the medium 
rest frame (q = 0). Since there seems to be quite a deal of confusion in the 
literature, we are are trying to be very explicit in this paper, showing the 
calculation of meson correlators in a Fermi gas of nucleons at q = 0. 

So far, the majority of the applications of the QCD-sum-rule method has 
dealt with the vector mesons. In this paper, having in mind possible future 
applications to other mesons, we collect the formulas for different channels as 
well. 



2 Dispersion relations and QCD sum rules 

QCD sum rules [26,27] for meson correlators at finite temperature [14-21] and 
finite density [22-25] have provided useful information on medium modifica- 
tions of mesonic spectra. The basic quantity in these studies is the retarded 
in-medium correlator of two meson currents, J and J', defined as 

n(w, q) = i I d'xe^'^eit) ({[J(x), J'(0)]» , (1) 

where ((...)) denotes the Gibbs average. Since the medium breaks the Lorentz 
invariance, the correlator depends separately on the energy variable, u, and 
the momentum, q, measured in the rest frame of the medium. The correlator 
(1) satisfies the usual fixed-q dispersion relation [17,19], 

n(.,q) = i/V , f-"' (2) 
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with p(y, q) = ImII(z/, q) denoting the spectral density. The RHS is not well- 
defined and the relation requires subtractions in order for the integral to con- 
verge. In QCD sum rules one applies the Borel transform, L Mg [26,27], to both 
sides of Eq. (2), where M B denotes the Borel mass parameter. The procedure is 
carried as follows: first the dispersion relation (2) is continued to the Euclidean 
space by replacing u with in, and then the Borel transform is carried in the 
variable f2 2 . As a result the RHS of Eq. (2) becomes -^pr J dv 2 e~ u l M Bp{y, q), 
which is well-defined. Thus, in fact, one does not need to perform subtractions 
in Eq. (2), since the Borel transform makes the spectral integral convergent. 
However, one can still perform the subtractions and in fact it is advantageous 
to do so in order to obtain more reliable sum rules for vector mesons both in 
the vacuum [26,27] and in medium [25]. The reason is that subtractions reduce 
the contribution from the high-lying resonances and continuum states in the 
spectral integral, as can be seen from the formula given below. 

The existent calculations of in-medium vector mesons [18,22,23] have implic- 
itly made use of the following subtracted dispersion relation, which follows 
from Eq. (2) when the subtraction is made at the point cu — |q| 

n(o;,q) -n(|q|,q) = 1 J , p(u, q) 

cu 2 -q 2 TV J (z/ 2 -q 2 )(z/ 2 -Cc> 2 )' { ) 



(We disregard the infinitesimal imaginary is terms here, since they are irrele- 
vant in QCD-sum-rule applications.) Compared to the study in the vacuum, 
there are two important differences. The first one is the appearance of the 
n(|q|, q) term, which is zero in the vector meson correlator in the vacuum, but 
in general does not vanish in medium (for instance, this is the case of the trans- 
verse vector meson correlator described in Section 4). The second difference 
is related to the analytic structure of the spectral density p(u, q) in nuclear 
medium. In addition to the singularities in the time-like region [y > |q|), the 
presence of medium induces singularities in the space-like region, [y < |q|) 
[28]. For instance, the particle-hole excitations of the nucleon Fermi sea at 
vanishing temperature lead to a cut reaching in the region < v < \q\kp/Ep, 
where hp and Ep are the Fermi momentum and energy. The spectral density 
connected with this cut is denoted by p sc {y, q). This phenomenon is called the 
scattering correction, or the Landau damping contribution, and its relevance 
for QCD sum rules was first pointed out by Bochkarev and Shaposhnikov [14] 
for the case of p mesons propagating in a pion gas. It has been subsequently 
used in Refs. [15-18,21-24]. 

The contribution of the scattering term is subtle in the limit of q — > 0, which 
describes the mesonic excitations at rest with respect to the nuclear medium. 
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The point is that in this limit one finds 

q^o 1/2 _ q 2 



where i? is a constant depending on the channel. The following sections provide 
a detailed analysis of this point. Thus, in the limit of q — > the following 
Borel-improved sum rule results from the subtracted dispersion relation (3) 



J M 



' n(^,0) \ _ lim q ^ n(|q|,q) R J_ T 2 -^/M| pM) 



where we have explicitly separated out the scattering term, and introduced 
p(v, 0) = p{v, q) — p sc {v, q). Note that if (and only if) Il(|q|, q) vanishes, then 
Eq. (3) may be viewed as the unsubtracted dispersion relation for the function 
n(o;,q) = U(uj , q) / \oo 2 — q 2 ). This is the case of the vector correlator in the 
vacuum, or of the longitudinal vector correlator in medium. But this is not 
the case in other channels, in particular, IT(|q|,q) does not vanish for the 
transverse vector channel in medium. In the following sections the issue will 
be discussed in detail on an example where the medium consists of the Fermi 
gas of nucleons, and the meson correlators are evaluated at the lp-lh level. 
This case is relevant, since existent QCD-sum-rule calculations at finite density 
use this approximation to model the scattering term in the phenomenological 
spectral function. A similar analysis can be carried for more complicated many- 
body treatment, and for other contributions to the correlators, such as the 
meson loops. 



3 Meson correlators in the Fermi gas of nucleons 



From now on our medium is the nuclear matter described by the Fermi gas 
of nucleons. Using the imaginary-time formalism [29,30] we evaluate the in- 
medium part of the meson correlator at the lp-lh level. A straightforward 
calculation leads to the formula 



U ab (q) = -J 



d 3 k 



2£ fc (27r) 3 
1 



e P{k-u-n) _|_ I 



Tab(k,q) + T ab (k- q,q) 



q 2 + 2k-q q 2 - 2k ■ q 



e f3(k-u+n) _|_ l 



T ab {-k -q,q) T ab (-k, q) 



q 2 + 2k-q 



2k ■ q 



(6) 



k°=E k 



Here q is the meson four-momentum, k is the internal four-momentum, k = 
E k = \JM 2 + k 2 , M is the nucleon mass, is the four- velocity of the medium, 
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Table 1 



List of channels considereded, with the corresponding vertices T a ,Ti„ and projection 
tensors V?t 



(3 is the inverse temperature, /x is the baryon chemical potential, and a and b 
label the channel. The Dirac trace factor is given by 



T ab (k, q) = tr \(f + M)T a {^ + j + M)V t 



(7) 



where T a and T b denote the appropriate vertices (cf. Table 1). Note that in the 
rest frame of the medium q = (v,q), where v is a purely imaginary bosonic 
Matsubara frequency, v = 2irin/[3, hence the denominators in (6) are well 
defined. At a later stage of the calculation we shall perform the standard 
analytic continuation in order to obtain the imaginary part of the retarded 
correlator. This is achieved by the substitution v = 2nin/ j3 — > vr + ie, where 
v R is real. 

Except for the scalar and pseudoscalar channels we need to project expression 
(6) on the appropriate quantum numbers using projection tensors Vfy: 

<pba 

n ab = E n (0 Pg, n (i) = n a6 ^-|-, di m p w = £ vfo c . (8) 



5 



The considered choices for T a , T b and the corresponding Vfy are listed in 
Table 1. We analyze all diagonal channels and one non-diagonal case: the 
vector-tensor channel. The explicit expressions for V?k are given in the Ap- 
pendix. The scalar functions V(i) depend on two Lorentz scalars: q 2 = v 2 — q 2 , 
and q ■ u = v. 

In the considered cases the trace factor (7) satisfies the following symmetry 
relations^: 



T ab (k - q,q) =T ab (k,-q) , T ab (-k - q,q) = T ab (k,q), 

T ab (-k,q) = T ab (k,-q). (9) 

This feature allows us to rewrite Eq. (6) in the following form 

** = f (10) 



where 



and 



= { 2k ■ q [T ab (k, q) - T ab (k, -q)] 

T- 

dimP(j) 



<pab 

- q 2 [T ab (h, q) + T ab (k, -<?)]} ^-S-- (12) 



We can decompose 



Afu) = Au) + B®x + C {i) x 2 , (13) 



where x is the the cosine of the angle between the three- vectors k and q, and 
the functions Ar^, Brq and Cu\ depend only on |k|, v and |q| . For simplicity of 
notation we shall from now on denote |k| and |q| by k andq, respectively. The 
explicit form of the coefficients A^),B^ and C(j) for the considered channels 
is given in Table 2 . 

The angular integration in (10) is elementary and yields 

1 This is not a universal feature and it does not hold for certain non-diagonal 
channels, which we do not take into consideration in this paper (e.g. the scalar- 
vector channel). 
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Table 2 

Explicit forms of the coefficients A^,B^ and C(j) denned by Eqs. (12) and (13), 
here q 2 = q • q and k 2 = k • k. 



j,2 _ ? 2)2 _ 4(_E feZ/ _ fc gx )2 



+4) ^ 



2£ fc z/ + z/ 2 - g 2 - 2fcg 



2£ fc z/ + z/ 2 - q 2 + 2kq 



+ In 



■BC(i)kq 

2E k v -v 2 + q 2 + 2kq 
2E k v - v 2 + q 2 - 2kq 



with 



± C (l) \2E k v ±u 2 T q 2 } 2 + 2kq [2A {l) kq + B (l) (2E k v ± v 2 T q 2 ) 



(14) 



(15) 



By substituting (14) in (10) we express IT^ as a single integral over k. 



4 Subtraction constant U(q, q) 



Equipped with explicit formulas from the preceding section we may now eval- 
uate the subtraction constant U(q,q). By inspection of Table 2, we find that 
the ratio J\f^(k, v = q, q)/(E k q — k • q) 2 is a number depending only on the 
channel. It vanishes in vector and axial-vector longitudinal channels, and in 
transverse tensor channels. In other cases Il(g, q) is given as an integral over 
the distribution function (11). For example, in the transverse vector channel 
we have 

1 r Hkk 2 

U VT (q,q) = - ^f(E k ,^,P). (16) 

7T J E k 

At zero temperature the integral over the distribution function can be done 
analytically with the result 



kp 

n%,g) = -i/-< 



d 3 k Nu)(k, v = q, q) 



2E k (2nf (E k q-k-qy 



7T" 



ir 2 E F 



where 



$(v F ) =v F + \\n 



i, 1-vf 

I-I+Vf 



(17) 



(18) 



and vf = kp/Ep = kp/^jM 2 + kp is the velocity of nucleons on the Fermi 
surface. Note that in our calculation IT^ (q, q) does not depend on q. The 
coefficients characterize the channel. They are given in Table 3. At low 
density $ can be expanded in powers of kp 

i k 3 3A- 5 

*M = Hkp/^P + ¥ F ) ~ + — ^, (19) 



which we will use explicitly in the following sections. 



5 Spectral density in the long-wavelength limit 



The spectral density acquires contributions in the time-like region (q > u), 
and also in the space-like region (q < v). The latter region is strictly related 
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to the presence of the medium. Below we shall not be concerned with the 
modifications of the time-like production cut by the medium. Since this cut is 
far away (in the Fermi gas and for q = it starts at v — 2E F ) its contribution 
to Eq. (5) is small. In addition, the limit of q — > is regular in the time-like 
region. This is not the case for the space-like region, where special care is 
needed when q — > [14]. The finite contributions stemming from the space- 
like region are known as the scattering terms, or Landau-damping terms. In 
this section we analyze in detail the q — > limit. Our results for scattering 
terms for various channels are shown in Table 3. 

Since we have adopted the imaginary-time formalism in our approach, we 
have been implicitly dealing with purely imaginary energies so far. Now, we 
are going to do the analytic continuation to real energies. This can be done 
in the standard way, by making the replacement v — > vr + is. In this way one 
obtains the retarded correlator for real energies ur. The imaginary part of the 
retarded correlator is generated by the two logarithms appearing in Eq. (14). 
The imaginary part of the first logarithm is 



7T 



9(-E k + y/E* + 2kq + q 2 - v) - 9{-E k + s[e 2 - 2kq + q 2 - u) 



+6(-E k - J E 2 k + 2kq + q 2 -u)- 9(-E k - J E 2 - 2kq + q 2 - v) 



(20) 



where we have omitted the index R indicating that v is real from now on. For 
space-like momenta in the long- wavelength limit {y < q and M » u, q — > 0) 
only the first two terms contribute and the imaginary part can be reduced to 
the expression 



v \ AM 2 a 
n0(k-\k4v,q)\), fc>,9)=2V?T^ + 1 -§- (21) 



In a similar way we deal with the second logarithm in (14) and find that its 
imaginary part in the discussed region is 



v AM 2 a 
n9(k-\k + (v,q)\), k+(v,q) = -^^—^ + 1 + 1. (22) 
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Consequently, the imaginary part of H^> is equal tcf 2 ] 



imn (i) 



1 



128vrg 3 



=-\ \k+\ 



dk 



(23) 



For space-like momenta in the long-wavelength limit both v and q tend to zero 
but with their ratio is kept constant 



q, v — > 0, u/q = a (0 < a < 1). 



(24) 



In this case we obtain 



aM q 



a 



- r 2 2 



(25) 



and the leading term of Eq. (23) is 



Imn (i) (a) = lim 



q^o 128vrg 3 



dk 
~EM' 



f(E k , fi, (3) C^{k, aq, q) - C (i) {k, aq, q) 



+ 



Q 



2E k k 



f(E K ,n,fl) \C^ } (K } aq,q) + C (l) (K,aq,q) 



(26) 



Of special interest is the case of cold matter, where (3 —>■ oo. In this situation 
the distribution function (11) reduces to the step function 9(fi — Ek) = 9(kp — 
k), and formula (23) becomes an integral over finite range of k 



ImII w (^<?) 



1287rg 3 



dk 



kp 



\k+\ 



-9(k F -\k+\) J — £^(A;,z/,g) 



(27) 



2 It is still possible to change the variables in (23) in such a way that the lower limit 
of integration in the two integrals is the same (E^ = ^(qx — v) and E\. = ^(qx + v) 
in the first and second integral, respectively.) This leads to the compact integral 
representation of the imaginary part valid for arbitrary temperature and chemical 
potential. For vector correlators using the appropriate forms of C v we recover the 
result of Bochkarev and Shaposhnikov [14] in this way. 
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In all the cases we consider turns out to be a simple polynomial and 
the integrals over k in (27) can be easily performed. We do not present the 
results of such integrations here, since the final expressions are rather lengthy. 
Moreover, we are concentrated on the long-wavelength limit, where further 
simplifications are possible. In this case 

6(k F -\k T \) = e(v F -^) = 9(v F -a), (28) 

and the zero-temperature limit is given by Eq. (26) with /3) replaced 

by 9(v F — ol). 

We note that lmH^\u, q)/(v 2 — q 2 ) becomes proportional to the delta function, 
S(u 2 ), in the limit (24), in agreement with Eq. (4) [14]. This can be seen 
explicitly from the integral 

h/^Sf =/|^>n<-( a )/(0) s fl ( „/(0), (29) 

where f{v 2 ) is an arbitrary continuous function and describes the strength 
of the distribution. Our results for various channels obtained for the zero- 
temperature case are shown in Table 3. 



6 Results for the vector channel 

In this section we discuss in some greater detail the case of the vector channel. 
This channel is particularly important in view of the experimental evidence of 
in-medium modification of light vector mesons. First, we recall the well-known 
fact that on general grounds at q = the transverse and longitudinal vector 
correlators are equal [28] 

nv T M) = n Vj >,o). (30) 

In order for this to hold, we note from Eq. (5) that we must have 

lim (IM^) -U Vl (q,q)) = - ^ . (31) 

Below we verify explicitly that indeed this equality holds for the Fermi gas at 
finite temperature. 
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The coefficients Au\, Bu\ and Cu\ for the vector channel are given in Table 2. 
After inserting their explicit form into Eq. (15) we find 

£$ L = -16k 2 (u 2 - q 2 ){AE 2 + u 2 -q 2 ± AE k u) (32) 



and 



1~Vt 



-lQk 2 (2k 2 q 2 - 2E 2 k v 2 T 2E k v{y 2 - q 2 ) - V — + O . (33) 



Substituting these equalities into Eq. (26) and (29) yields, after elementary 
algebra, the following expression: 



Rv T — Rvt 3 



7T 





— J a 2 da J kdkf(E k ,n,f3) 



+ ^Jj^y^ m ^' P) ' (34) 



where k = aM/yl — a 2 , and E K = M/yl — a 2 . Next, we interchange the 
order of the a and k integrations in the first term in Eq. (34), and replace 
the integration variable from a to k in the second term. As a result we obtain 
exactly the integral for lim^o Hy T (q, q), as given in Eq. (16). Since Uv L (q, q) = 
0, it completes the explicit check of Eq. (31) in the Fermi gas. 

Let us now turn to the discussion of the zero-temperature case. Through the 
use of Eq. (32) in (27) we recover Eq. (5.75) of Ref. [28] 

ImlT^ (a) = ~ <* 2 W"F ~ <*) E 2 (35) 

2ir 



The strength of the distribution function as defined in (4) or (29) is 

= "si;' < 36 > 



On the other hand, the explicit calculation in the transverse vector channel 
gives 

* 7 . . a(a 2 — v 2 p)6(vf — ot) ^ 9 , „. 

Imn yT (a) = -i -E F (37) 

An 
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and the corresponding strength is 



R 



V T 



1 

6n 



hip 



(38) 



The calculation of the imaginary part in other channels can be done in the 
analogous way. In Table 3 we summarize the results for all the channels giv- 
ing the exact and approximate strength I&\ and the constant which 
determines the subtraction term U.^(q, q). 

To conclude this section, we comment on the use of the scattering term in 
the existent calculations concerning vector mesons in medium. For symmetric 
nuclear matter hp = 3ir 2 p/2, where p is the baryon density. Hence, in the limit 
(24) one can write — |lmIly L /(^ 2 — q 2 ) — > p sc 5{y 2 ) / (Sit), where p sc = 2it 2 p/E F . 
The last result agrees with the normalization of the scattering term used in 
the literature for both the longitudinal and transverse channels [16-18,21-24]. 
As we have shown, the scattering term for the transverse and longitudinal 
spectral densities is different. However, the presence of the subtraction con- 
stant lim^o ny T (q, q) compensates this difference and, in effect, the "cavalier" 
approach of Ref. [16-18,21-24] indeed leads to correct results. 



7 Summary 



In this paper we have calculated the scattering term contributions, as well as 
corresponding subtraction constants in various mesonic channels. We stress 
that these subtraction constants are well-defined (unambiguous) in a given 
many-body treatment. This is in contrast to the case of the vacuum, where 
the subtraction constants, connected to renormalization, are not well-defined. 
However, the presence of the medium does not bring new divergences to the 
theory; in other words, medium effects are finite. This is reflected by the 
finiteness of our constants U(q,q). We have also shown the necessity of the 
subtraction constants for the consistency of the approach: together with the 
scattering terms they lead to the required equality of the correlators in the 
transverse and longitudinal channels when q — > 0. 

As is well known, for finite q the longitudinal and transverse correlators are no 
longer equal. Our expressions for the scattering term can be straightforwardly 
extended to that case. The subtraction constants do not depend on q in the 
Fermi-gas model. 

The explicit calculations of this paper were done for the Fermi gas of nucleons. 
In Table 3 we show our results, needed for the construction QCD sum rules in 
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Table 3 



List of our results. The scattering terms are characterized by the strength R defined 
by Eq. (29), whereas the subtraction constants are characterized by the coefficient 
W defined in Eq. (17). The function <l> is introduced in Eq. (18), and the limit of R 
at low density is obtained from Eq. (19). 

these channels. We note, that the calculation in the Fermi gas is a model calcu- 
lation. Other effects, such as meson production, rescattering, meson-exchange, 
etc. will result in modifications of both the scattering terms and the subtrac- 
tion constants. These effects need not be small and should be incorporated. 
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Tensors L^, and are defined as follows 3 



t- = <r - - w-v™"W-v™ v ) i (39 ) 

q ■ q — q ■ u 2 



3 Note that our definition of the T and L tensor has the sign chosen is such a way 
that they are projection operators. Opposite convention is frequently used, in order 
to make ImIly L /(^ 2 — q 2 ) and Imnv T positive. 
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q ■ q q ■ q — q ■ u 2 

Q^ = ^C (41) 

q-q 

The two projection operators L^ v and T Miy are both transverse to q. In addition, 
rpnu j g transverse to u^, while is longitudinal to u**. The projection operator 
Q^ u is longitudinal to g. It appears in the decomposition of the axial vector 
correlator, since the axial-vector current is not conserved. Tensors L^ u and T^ v 
can be obtained as a sum over polarization vectors [31]. Eqs. (39-41) imply 
the following projection-operator relations: 

t«tj = t;, l«l: = l;, q«q: = q», (42) 

t«l: = t«q: = q«l: = o, (43) 

T/ = 2, L; = l, Q; = i. (44) 



The covariant decomposition of the vector and axial-vector correlators into 
tensors (39) - (41) follows directly from the Lorentz structure of expression 
(6). In the rest-frame of the medium, = (1,0,0,0), our definitions of the 
projection tensors reduce to commonly used non-covariant expressions. 

We also give the appropriate projection tensors for other channels, involving 
the tensor coupling. Such correlators have been considered in the vacuum 
[32,33]. Our expressions prepare grounds for the extension of such calculations 
to finite density and temperature. The tensor correlators can be decomposed 
in the basis of the projection tensors T^ a/3 and L^] a/3 , where the sign (±) 

refers to the parity of the excitation. The tensors L^~J. a/3 and L^). a ^ have the 
following structure 



L iJ; a p = ^2 [V* W + L ^q^» - L ^quq a - L va q^q p \ , (45) 
= ^ e nvar£ a pK\L TK q a q x , (46) 

where is defined in (40). Tensors T^ a/3 and T^ aj3 follow from equations 
(45) and (46) with L^ u replaced by T^ v . In analogy to Eqs. (42-44) we find 

rp(±) rpaP\(TT _ rp(±) UT j (±) j a/3;ar _ j (±) err ( aij\ 

1 ^u;al3 1 (±) — 1 i iu; > %i/;a/3^(±) M^; ' ) 

TW^ = 2, Lg^ = l. (49) 
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We note that all products of the tensors with opposite parity vanish. 

Tensors L M; a/3 and T^. a/3 appear in the decomposition of the correlators in the 
vector-tensor channel. They are defined as 



L^, a /3 = y ^2 (V*<7/3 ~ L »pQa) , (50) 



and 



T tx\ a/3 — y ^2 (^"9/9 ~ ^/3?a) • (51) 



The normalization of L^. Q/ g and T^. a p is induced by the requirement that the 
algebra of projection tensors is closed. For example, the products of L u . a/3 with 
L l and L CT l fi are 

L^apL^' UT = I/ m;ctt , L^L v . a p = L^ap. (52) 
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